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ABSTRACT 

A knowledge  of  the  brightness  of  moonlight  is  needed  for  detailed  calculations  of  the  limiting 
magnitude  of  astronomical  detectors,  whether  they  be  visual,  photographic,  or  electronic.  The 
previous  literature  contains  no  method  for  making  even  crude  estimates  and  has  few  actual 
measurements  of  moonlight  brightness.  In  this  paper  we  report  new  measurements  of  the  sky 
brightness  from  the  2800-m  level  of  Mauna  Kea.  In  addition,  we  present  a model  for  predicting  the 
moonlight  as  a function  of  the  Moon’s  phase,  the  zenith  distance  of  the  Moon,  the  zenith  distance  of 
the  sky  position,  the  angular  separation  of  the  Moon  and  sky  position,  and  the  local  extinction 
coefficient.  The  model  equations  can  be  quickly  calculated  on  a pocket  calculator.  A comparison  of 
our  model  with  our  lunar  data  and  with  some  Russian  solar  data  shows  the  accuracy  of  our 
predictions  to  range  from  8%  to  23%. 

Key  words:  night-sky  brightness-moonlight-atmospheric  extinction 


1.  Introduction 

It  is  of  fundamental  interest  to  know  the  limiting  mag- 
nitude of  visibility  for  finding  objects  with  a telescope. 
This  subject  has  recently  been  addressed  by  Schaefer 
1990a  based  on  a database  of  314  visual  observations  of 
stars  viewed  in  small-  to-medium-sized  telescopes  with  a 
variety  of  magnifications.  For  large  telescopes  using  video 
acquisition  and  guiding  systems,  the  question  is  still  im- 
portant. It  is  also  relevant  to  the  prediction  of  proper 
exposure  times  for  optical  CCD  cameras.  A related  ques- 
tion is:  what  is  the  magnitude  limit  when  the  Moon  is  out? 
With  the  eye,  video  system,  or  CCD  camera  the  detec- 
tion of  an  object  depends  on  the  contrast  of  the  brightness 
of  the  object  in  question  versus  the  brightness  of  the  sky. 
We  can  think  of  this  as  achieving  a certain  signal-to-noise 
ratio.  Whatever  brightens  the  sky  degrades  the  limiting 
magnitude  of  the  eye  or  video  system  by  approximately 
the  same  amount. 

In  a recent  paper  (Krisciunas  1990)  we  presented  a 
sample  of  photometric  measurements  of  the  sky  bright- 
ness, with  and  without  moonlight.  Most  of  the  measure- 
ments (and  all  of  the  moonlight  observations)  were  ob- 
tained at  the  2800-m  level  of  Mauna  Kea.  The  telescope  is 
a 15-cm  Newtonian  reflector,  and  the  photometer  em- 
ploys standard  UBV  filters  and  an  uncooled  RCA  931 A 
photomultiplier.  The  beam  size  is  6.5  square  arc  minutes. 

At  a site  not  hampered  by  light  pollution  it  is  found  that 


the  inherent  dark  sky  level  at  the  zenith  (but  away  from 
the  galactic  plane  and  more  than  two  hours  after  the  end 
of  twilight)  is  about  V — 21.8  to  22.0  mag  sec-2,  B = 22.8 
to  23.0  mag  sec-2.  At  solar  maximum  these  inherent 
levels  are  0.4  to  0.6  mag  sec-2  brighter. 

The  lunar  sky  brightness  is  the  arithmetic  difference 
between  the  observed  sky  brightness  with  the  Moon 
above  the  horizon  and  the  inherent  dark  sky  value,  given 
the  phase  of  the  solar  cycle.  In  Figure  8 of  Krisciunas 
1990  these  differences  are  plotted  for  21  measurements, 
where  the  nominal  dark  sky  levels  are  taken  to  be  the 
yearly  averages  for  the  zenith.  A brief  summary  of  some 
measurements  is  also  given  by  Walker  1987  based  on 
measurements  at  Cerro  Tololo.  Walker’s  measurements 
made  within  4 days  of  full  Moon  were  taken  90  degrees 
or  more  away  from  the  Moon. 

The  effect  of  the  Moon  on  the  sky  brightness  at  some 
position  is  a function  of  the  phase  of  the  Moon,  the  zenith 
distance  of  the  Moon,  the  zenith  distance  of  the  sky 
position,  the  angle  between  the  Moon  and  the  sky  posi- 
tion, and  the  atmospheric  extinction.  In  this  paper  we 
report  additional  moonlight  brightness  observations  from 
the  2800-m  level  of  Mauna  Kea.  In  addition,  we  derive  a 
simple  model  that  gives  the  relation  between  the  scat- 
tered moonlight  and  the  five  variables  mentioned  above. 
This  model  results  in  an  easy-to-use  equation  that  will 
predict  the  moonlight  brightness  with  an  accuracy  of 
8%  to  23%. 
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2.  Observations 

A dozen  observations  of  the  brightness  of  the  moonlit 
sky  were  made  from  the  2800-m  level  of  Mauna  Kea  with 
the  identical  instrumental  setup  as  used  by  Krisciunas 
1990.  These  dozen  observations  are  to  be  joined  with  the 
21  observations  in  Table  5 of  Krisciunas  1990.  In  Table  1 
we  give  the  new  measurements  of  the  observed  V-band 
sky  brightness  on  moonlit  nights.  The  first  six  columns 
give  the  UT  date  and  UT  time  of  the  observation,  the 
azimuth  and  zenith  distance  of  the  sky  position,  and  the 
azimuth  and  zenith  distance  of  the  Moon.  Note  that  all 
angles  have  been  rounded  to  the  nearest  degree  since 
greater  accuracy  is  not  needed.  The  seventh  column  gives 
p,  which  is  the  angular  distance  between  the  Moon  and 
the  sky  position.  The  eighth  column  gives  the  phase  angle 
(a)  of  the  Moon,  which  is  defined  as  the  angular  distance 
between  the  Earth  and  the  Sun  as  seen  from  the  Moon. 
(When  the  Moon  is  between  full  and  new,  a < 0.)  The 
phase  angles  were  calculated  using  the  simpler  algorithm 
of  Meeus  1980  (eq.  (31.4)),  which  ignores  the  noncon- 
stancy of  the  (ecliptic)  latitude  of  the  Moon.  The  ninth 
column  gives  the  observed  V magnitude  of  the  sky  as 
measured  in  magnitudes  per  square  arc  second. 

The  last  column  in  Table  1 gives  the  sky  brightness  of 
the  moonlight  alone  as  measured  in  nanoLamberts  (nL). 
The  reason  for  converting  from  a logarithmic  unit  (magni- 
tudes per  square  arc  second)  to  a linear  unit  (nanoLam- 
berts) is  because  a magnitude  difference  compares  the 
moonlight  to  the  background  which  is  variable  from  site  to 
site.  In  addition,  the  use  of  a linear  unit  allows  for  the  easy 
subtraction  of  the  background  sky  brightness.  The  choice 
of  the  nanoLambert  as  the  linear  unit  is  to  follow  the  lead 
of  Knoll,  Tousey  & Hulburt  1946,  Weaver  1947,  Hecht 
1947,  Garstang  1989,  and  Schaefer  1990a.  The  relation 

TABLE  1 


V-band  sky  brightness  on  nights  with  moonlight3 


UT 

Skv  DC 

sn. 

Moon 

p 

a 

V 

Bmoon 

(Obs.) 

AZ 

Z 

AZ 

Zm 

Date 

Time 

O 

O 

O 

o 

O 

O 

mag/ sec 2 

nL 

1990 

12  Feb 

7:11 

280 

9 

94 

79 

88 

-30 

19.813 

308 

28  Dec 

7:05 

89 

18 

300 

4 

22 

49 

18.907 

833 

1991 

4 Jan 

8:07 

125 

17 

85 

82 

68 

-48 

20.141 

206 

4 Jan 

8:15 

63 

38 

85 

80 

45 

-48 

19.605 

380 

4 Jan 

8:23 

55 

30 

86 

78 

53 

-48 

19.684 

354 

4 Jan 

8:30 

82 

52 

87 

77 

24 

-49 

18.615 

1087 

4 Jan 

8:34 

143 

13 

87 

76 

69 

-49 

19.902 

281 

4 Jan 

8:42 

99 

61 

87 

74 

16 

-49 

18.456 

1255 

4 Jan 

8:49 

79 

68 

88 

72 

9 

-49 

18.228 

1557 

4 Jan 

8:54 

165 

11 

88 

71 

69 

-49 

19.771 

329 

21  Jan 

6:26 

263 

49 

265 

62 

14 

119 

20.572 

71 

21  Jan 

6:33 

113 

24 

266 

64 

85 

119 

21.031 

30 

Observations  made  at  2800-m  elevation  of  Mauna  Kea,  Hawaii. 


between  the  sky  brightness  in  nanoLamberts  ( B ) and  in 
magnitudes  per  square  arc  second  (V)  is 

B = 34.08  exp(20.7233  - 0.92104  V)  , (1) 

as  given  in  equation  (27)  of  Garstang  1989.  The  surface- 
brightness  equivalent  to  one  tenth  magnitude  star  per 
square  degree  (that  is,  the  S10  unit)  equals  0.263  nL.  An 
alternative  expression,  equivalent  to  equation  (1),  which 
includes  the  definition  of  the  S10  unit  and  the  scale  factor 
to  convert  to  nanoLamberts,  is 

B = 0.263  a(<?-v)  , 

where  a = (100)02  « 2.51189  and  Q = 10.0  + 2.5 
log(36002)  ~ 27.78151. 

The  dark  nighttime  sky  brightness  (B0)  as  a function  of 
zenith  distance  (Z)  is 

Bo(Z)  = Bzen10-°'4W-1)X  , (2) 

where 

X = (1  — 0.96  sin2Z)-0'5  . (3) 

The  extinction  coefficient  is  k in  units  of  magnitudes  per 
air  mass.  X is  the  optical  pathlength  along  a line  of  sight  in 
units  of  air  masses.  Note  that  the  formula  for  X in  equation 
(3)  is  appropriate  for  the  night  glow.  Equations  (2)  and  (3) 
are  simplifications  of  equations  (29)  and  (30)  from 
Garstang  1989  that  are  adequate  for  the  present  purposes. 
The  dark  zenith  sky  brightness  (Bzen)  is  taken  as  the  yearly 
average  reported  in  Krisciunas  1990  or  as  discussed  be- 
low. The  extinction  coefficients  were  taken  to  equal  the 
value  measured  for  that  night  or  the  median  value  at  the 
2800-m  level  of  Mauna  Kea,  namely  0.172  mag/air  mass, 
obtained  by  Krisciunas  1990. 

The  last  column  of  Table  1 has  the  background  sky 
brightnesses  subtracted  out.  That  is,  the  value  for  the 
background  alone,  B0  (Z),  is  subtracted  from  the  observed 
brightness  for  the  moonlight  plus  background,  B , leaving 
the  brightness  of  the  moonlight  alone,  Bmoon.  When  the 
moonlight  is  dim  compared  to  the  background,  the  natu- 
ral variations  in  the  background  will  cause  large  relative 
errors  in  derived  values  of  Bmoon.  Therefore,  we  have 
excluded  nine  observations  (for  which  the  moonlight  is 
fainter  than  the  dark  sky  background)  from  our  statistical 
analysis  of  the  model  compared  to  actual  observations. 

The  observations  of  1991  January  4 (UT)  are  particu- 
larly useful,  since  a concerted  effort  was  undertaken  to 
narrow  down  the  functionality  of  the  components  of  the 
lunar  sky-brightness  effect.  This  involved  obtaining  a 
dark-sky  value  on  that  night  (V  = 21.444  mag  sec-2)  an 
hour  before  the  Moon  rose.  This  value  compares  well 
with  the  1990  average  (the  year  of  sunspot  maximum)  of 
V = 21.377  ± 0.095  mag  sec-2,  based  on  measurements 
of  5 nights.  For  the  two  observations  in  Table  1 from  1990 
we  shall  adopt  the  1990  yearly  average  as  the  nominal 
dark-sky  value.  We  shall  adopt  the  actual  dark-sky  value 
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for  1991  January  4 for  reference  in  calculating  the  lunar 
effect  on  that  night.  On  this  night  we  also  measured  the 
V-band  extinction  to  be  k = 0. 116  ± 0.019  mag/air  mass, 
based  on  8 measurements  of  5 stars  over  a range  of  1.5  air 
masses.  For  the  1991  January  21  data  we  shall  adopt  a 
dark-sky  value  of  V = 21.388  mag  sec-2,  equal  to  the 
average  of  the  six  measurements  of  the  previous  year. 

3.  Atmospheric  Scattering 

The  atmospheric  scattering  of  moonlight  is  a complex 
problem.  The  sky  brightness  will  depend  on  the  extinc- 
tion coefficient,  the  magnitude  of  the  Moon,  the  zenith 
distance  of  the  Moon,  the  angular  separation  between  the 
Moon  and  the  sky  position,  and  the  zenith  distance  of  the 
sky  position.  In  these  next  two  sections  we  describe  a 
simple  scattering  model  (elaborated  by  one  of  us,  B.  E.  S.) 
so  as  to  establish  functional  forms  for  the  variation  of  each 
parameter  that  are  fairly  accurate. 

When  looking  out  along  a line  of  sight,  an  observer  will 
detect  scattered  light  from  the  Moon.  Let  the  field  of  view 
be  a circle  with  an  angular  radius  of  c|),  so  that  the  scat- 
tered light  must  come  from  inside  a narrow  cone.  The 
volume  of  the  cone  can  be  broken  up  into  differential 
volumes 

dV='n(R<\>fdR  , (4) 

where  R is  the  distance  from  the  observer  to  the  volume. 
This  volume  will  have  some  density  of  scatterers,  D, 
which  will  be  proportional  to  the  volume  extinction  coeffi- 
cient b 

D = cb  , (5) 

where  c is  some  constant.  The  number  of  scatterers  in  the 
differential  volume,  dN , will  be 

dN  = D dV  = it  cb  R2$2dR  . (6) 

The  brightness  of  scattered  light  from  each  volume  will  be 
proportional  to  dN. 

The  differential  volume  will  be  illuminated  by  moon- 
light dimmed  by  the  atmosphere.  If  I*  is  the  illuminance 
of  the  Moon  outside  the  atmosphere  and  I is  the  illumi- 
nance of  the  Moon  inside  the  atmosphere,  then  these  two 
quantities  will  be  related  as 

1 = 1*  10-0-4*Xm  , (7) 

where  k is  the  extinction  coefficient  and  Xm  is  the  optical 
pathlength  to  the  Moon  as  measured  in  air  masses  (see 
below).  The  illuminance  of  the  Moon  (in  footcandles)  can 
be  related  to  the  V magnitude  of  the  Moon  (m)  as 

j*  _ 2Q_0-4(m  + 16-57)  (Qj 

(see  eq.  (16)  of  Schaefer  1990a).  The  magnitude  of  the 
Moon  will  be  a function  of  the  phase  angle,  a in  degrees, 
such  that 


m = -12.73  + 0.026  lal  + 4 X 10-9a4  (9) 

(see  page  144  of  Allen  1976).  The  brightness  of  scattered 
moonlight  will  be  proportional  to  I.  Equation  (9)  gives  m 
= — 12.73  for  full  Moon,  in  agreement  with  the  results  of 
Lane  & Irvine  1973  extrapolated  to  a = 0.  However,  this 
ignores  the  “opposition  effect”.  For  lal  < 7°  the  bright- 
ness of  the  Moon  deviates  from  this  relation  (see,  e.g., 
Whitaker  1969).  When  the  Moon  is  exactly  full  it  is  about 
35%  brighter  than  the  extrapolation  would  predict,  as- 
suming, of  course,  that  it  is  not  undergoing  a penumbral 
or  umbral  eclipse. 

The  luminous  intensity  of  the  differential  volume,  dL , 
will  be  proportional  to  the  number  of  scatterers  in  the 
volume,  the  illuminance  of  the  Moon,  and  the  scattering 
function.  The  scattering  function,  / (p),  describes  the 
intensity  of  scattered  light  as  a function  of  the  scattering 
angle  p.  Therefore, 

dL  = Z/(p)  dN  . (10) 

In  this  equation  and  below,  the  various  constants  of  pro- 
portionality will  be  absorbed  into  the  scattering  function. 
The  light  traveling  from  the  differential  volume  to  the 
observer  will  be  attenuated  not  only  by  the  atmosphere 
but  also  by  the  inverse  square  law  for  the  intensity  of 
light.  The  perceived  surface  brightness  from  the  differen- 
tial volume,  dB , will  be  the  perceived  brightness  divided 
by  the  subtended  solid  angle,  so  that 

dB  = dLe  Vir^R2  = //( p)  b e'T  dR  . (11) 

Here,  t corresponds  to  the  optical  depth  between  the 
observer  and  the  volume. 

The  total  apparent  surface  brightness  can  be  found  by 
integrating  along  the  entire  line  of  sight.  For  the  simple 
model  that  the  atmosphere  is  uniform  up  to  some  height 
H,  then  the  optical  depth  will  be  simply  bR.  The  integral 
must  run  from  zero  to  the  distance  of  the  top  of  the 
atmosphere  (that  is  H sec(Z)).  Thus,  the  surface  bright- 
ness will  be 

B = If(p)[l-e~bHsec{z)]  . (12) 

For  such  an  atmospheric  model,  the  extinction  coefficient 
can  be  related  to  the  volume  extinction  coefficient  as 

10-°Ak  = e~bH  . (13) 

The  air  mass  is  often  evaluated  as  the  secant  of  the  zenith 
distance,  and  this  is  accurate  for  any  direction  not  near  the 
horizon.  So,  the  secant  of  the  zenith  distance  in  equation 
(12)  will  be  better  approximated  as  the  air  mass  along  the 
line  of  sight,  X(Z). 

Many  expressions  for  air  mass  that  are  valid  low  on  the 
horizon  have  been  proposed  (Hardie  1962;  Rozenberg 
1966;  Garstang  1989;  Schaefer  1989).  Ideally,  each  extinc- 
tion component  must  be  handled  separately  since  they  all 
have  different  height  distributions  (cf.  Schaefer  1990b). 
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From  these  correct  calculations,  we  have  found  that  equa- 
tion III.  3,  17  from  Rozenberg  1966  is  the  most  reasonable 
of  the  simple  formulae  for  zenith  distances  all  the  way  to 
the  horizon.  Therefore, 

X(Z)  = [cos(Z)  + 0.025  e- 11  cos(z)]~x  . (14) 

In  this  formula  the  air  mass  goes  as  sec(Z)  when  far  from 
the  horizon  but  is  limited  to  40  air  masses  on  the  horizon. 
This  is  the  best  formula  to  use  for  correcting  the  bright- 
ness of  an  object  directly  observed  (e.g. , the  Moon  itself). 
However,  for  scattered  light  from  the  Moon  we  find  that 
air  masses  based  on  equation  (3)  give  the  better  fit  of 
observed  lunar  brightness  vs.  model  brightness,  for  a 
very  large  range  of  Bmoon.  We  might  call  equation  (3)  the 
“scattering  air  mass”  and  equation  (14)  the  “extinction  air 
mass”,  and  it  is  clear  from  our  data  that  equation  (14)  leads 
to  a gross  underestimate  of  the  lunar  sky-brightness  effect 
when  the  Moon  is  low  on  the  horizon. 

From  equations  (7),  (12),  and  (13),  the  model  surface 
brightness  will  be 

Bmoon  =/( p)  r 10-°  ^ [1-10  °™]  . (15) 

The  evaluation  of/(p)  will  be  given  in  the  next  section. 

This  last  equation  embodies  a number  of  complexities, 
the  details  of  which  have  not  been  fully  elaborated.  The 
most  obvious  are: 

(1)  The  volume  extinction  coefficient  h will  be  a sum  of 
two  exponential  functions  of  height. 

(2)  The  extinction  in  equation  (7)  will  be  a function  of 
the  height  of  the  volume. 

(3)  The  effects  of  a spherical  Earth  and  refraction  be- 
come important  near  the  horizon.  Better  calculations  can 
and  have  been  made  (cf.  Rozenberg  1966  and  references 
therein).  However,  the  results  are  of  much  greater  com- 
plexity. 

(4)  The  scattered  light  has  a small  but  significant  frac- 
tion of  multiply  scattered  light,  whereas  the  use  of  the 
scattering  function  in  equation  (12)  implicitly  assumes 
single  scattering. 

Equation  (15)  has  also  ignored  several  effects,  some  of 
them  quite  minor: 

(1)  The  Moon  has  an  asymmetric  distribution  of  maria 
so  that  the  Moon  and,  hence,  the  scattered  moonlight 
should  be  slightly  brighter  before  full  Moon  than  after 
(Stebbins  & Brown  1907;  Russell  1916;  Rougier  1934). 
However,  this  effect  is  very  small  and  Lane  & Irvine  1973 
were  unable  to  detect  it. 

(2)  For  lal  < 7°  the  lunar  opposition  effect  should  be 
taken  into  account,  in  which  case  the  derived  value  of 
Bmoon  must  be  multiplied  by  a factor  in  between  1.00  and 
1.35.  This  factor  only  comes  into  play  within  a day  of  full 
Moon. 

(3)  Refraction  will  change  the  altitudes  and  separation 
of  the  Moon  and  the  sky  position  from  those  used  in  this 


model.  However,  even  on  the  horizon,  the  altitude  will 
only  shift  by  half  a degree  at  most  and  will  change  the 
predicted  brightness  by  only  small  amounts. 

(4)  The  position  of  the  Moon  as  used  in  this  paper  is 
that  of  the  Moon’s  center,  whereas  it  would  be  more 
accurate  to  use  the  position  of  the  light  centroid.  How- 
ever, under  extreme  conditions,  the  difference  in  posi- 
tion is  a quarter  of  a degree  at  most,  so  that  this  effect  is 
negligible. 

(5)  The  Moon  is  sometimes  closer  and,  hence,  larger 
and  brighter  than  at  other  times  for  the  same  phase  angle. 
The  apparent  diameter  of  the  Moon  will  vary  by  11.6%  on 
average  from  apogee  to  perigee,  so  that  the  Moon’s  mag- 
nitude, as  derived  by  equation  (9),  may  be  in  error  by  as 
much  as  ± 0. 12  magnitude. 

An  improved  model  would  account  for  all  of  the  above 
effects  and  more.  Yet  the  real  atmosphere  is  never  as  well 
behaved  as  in  any  model  ever  calculated.  In  particular, 
the  concentration  and  distribution  of  the  atmospheric 
aerosols  are  highly  variable  with  time  and  location  (cf. 
Fig.  56  of  Rozenberg  1966).  With  such  large  variations,  it 
is  useless  to  refine  a model  intended  for  general  applica- 
tion. 

The  real  justification  for  the  simplicity  of  the  model  in 
equation  (15)  is  that  it  reproduces  the  observations  with 
remarkable  accuracy.  We  will  show  below  that  the  model 
reproduces  the  Mauna  Kea  data  (from  scattered  moon- 
light) with  an  rms  uncertainty  of  23%  and  the  data  (from 
scattered  sunlight)  from  two  Russian  sites  (Pyaskovskaya- 
Fesenkova  1957)  to  an  accuracy  of  8%  to  11%.  This  accu- 
racy is  achieved  despite  the  hundredfold  variation  in  the 
observed  brightness  of  scattered  moonlight. 

4.  Scattering  Function 

Equation  (15)  gives  a simple  formula  for  the  moonlight 
brightness,  where  the  functional  forms  for  /*,  k,  Xm,  and 
X are  derived  from  theory.  However,  the  scattering  func- 
tion, /( p),  has  not  yet  been  evaluated.  The  scattering 
function,  /( p),  is  proportional  to  the  fraction  of  incident 
light  scattered  into  a unit  solid  angle  with  a scattering 
angle  p.  The  scattering  angle,  p,  will  be  equal  to  the 
angular  separation  between  the  Moon  and  the  sky  posi- 
tion for  single  scattering. 

The  scattering  function  will  be  composed  of  additive 
terms  associated  with  the  two  types  of  scattering  in  our 
atmosphere.  The  first  type  is  Rayleigh  scattering  from 
atmospheric  gases,  which  will  contribute /R(p).  The  sec- 
ond type  is  Mie  scattering  by  atmospheric  aerosols,  which 
will  contribute /M(p).  The  two  terms  will  add  (see  eq.  II.  2, 
38  of  Rozenberg  1966),  so  that 

/(p)=/r(p)+/m(p)  • (16) 

Remember  that  the  scattering  functions  have  absorbed 
constant  factors  relating  to  unit  conversions  and  normal- 
izations. 
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For  Rayleigh  scattering  from  atmospheric  gases,  the 
contribution  to  the  scattering  function  will  be 

/R(p)  = CR[1.06  + cos2(p)]  (17) 

(Rozenberg  1966).  CR  is  a proportionality  constant  yet  to 
be  determined. 

For  Mie  scattering  of  aerosols,  the  contribution  to  the 
scattering  function  will  be  /M(p).  Explicit  equations  for 
/M(p)  have  been  compiled  in  the  treatise  of  van  de  Hulst 
1957  for  many  particle  shapes,  sizes,  and  optical  proper- 
ties. Unfortunately,  these  idealized  cases  have  no  utility 
for  the  real  atmosphere  where  the  aerosol  size  distribu- 
tion is  broad,  and  the  particle  shape  and  composition  are 
widely  variable.  In  principle,  /M( p)  can  be  found  by  inte- 
grating the  single  particle  functions  over  the  size,  shape, 
and  optical  properties  of  the  aerosols,  but  in  practice  this 
knowledge  is  never  known  with  sufficient  accuracy  to 
justify  the  integration.  Henyey  & Greens tein  1941  offer  a 
scattering  function  that  is  merely  a convenient  functional 
form  with  no  basis  in  data.  Thus,  the  only  alternative  is 
that  of  empirical  measurement. 

Therefore,  in  practice,  we  will  use  the  observations 
from  the  previous  section  to  define  the  scattering  func- 
tion. This  can  be  done  because  Bmoon,  I *,  k , Xm,  and  X 
are  known  for  all  observations  from  Mauna  Kea,  so  that 
equation  (15)  can  be  solved  for/(p).  Figure  1 shows  a plot 
of  the  scattering  function  as  a function  of  the  Moon/sky 
separation. 

Pyaskovskaya-Fesenkova  1957  presents  smoothed  data 
for  the  daytime  sky  brightness  that  can  also  be  used  to 
evaluate  /( p).  The  reason  is  that  scattered  sunlight  is 
identical  in  everything  but  intensity  with  scattered  moon- 
light. The  difference  in  intensity  can  be  easily  accounted 


Fig.  1-The  scattering  function,  /( p),  as  deduced  from  the  moonlight 
observations  reported  in  this  paper  and  in  Krisciunas  1990.  That  is,  for 
the  observations,  all  the  quantities  in  equation  (15)  are  known  so  that  the 
scattering  may  be  solved  for.  The  scattering  function  is  a function  of  the 
scattering  angle,  p.  The  scattering  function  from  equation  (21)  is  drawn 
as  the  smooth  curve. 


for  by  using  the  Sun’s  magnitude  (msun  = —26.74)  instead 
of  the  Moon’s  magnitude  in  equation  (8).  The  Russian  data 
are  for  the  whole  sky  from  sites  with  k = 0. 15  and  k = 0. 24 
and  for  solar  zenith  distances  of 0,  30,  60,  and  80  degrees. 
We  have  extracted  many  representative  points,  calcu- 
lated /(p)  from  equation  (15),  and  plotted  the  results  in 
Figure  2. 

The  first  item  to  notice  is  that  the  scattering  function  is 
identical  (to  within  the  scatter)  for  the  three  sites  with 
extinction  coefficients  of  0.15,  0.172,  and  0.24.  Second, 
there  is  no  systematic  difference  between  solar  and  lunar 
data.  The  third  and  fourth  items  are  that  the  data  with  the 
zenith  distance  of  the  source  and  of  the  sky  position 
greater  than  80  degrees  are  not  systematically  different 
from  the  average.  These  four  items  show  that  the  func- 
tional dependence  on  k,  I , Xm,  and  X is  correct.  This  is 
the  empirical  justification  that  the  approximations  leading 
to  equation  (15)  are  good. 

The  Mie  scattering  by  aerosols  is  highly  forward  scat- 
tering. Thus,  the  Rayleigh  scattering  will  dominate  for 
scattering  angles  greater  than  90  degrees.  Therefore,  the 
empirical  measures  of/(p)  for  large  angles  can  be  used  to 
measure  CR.  The  value  will  be  the  average  (over  data 
values  with  a scattering  angle  of  greater  than  90  degrees) 
of  the  empirical/(p)  divided  by  1.06  + cos2(p).  We  find  a 
value  of  2. 27  X 105forCR. 

The  Rayleigh  component  of/(p)  can  then  be  subtracted 
out,  leaving  only/M(p).  For  scattering  angles  from  10  to  80 
degrees,  the  empirical /M(p)  is  roughly  a straight  line  on  a 


0 30  60  90  120  150 

Scattering  Angle  (degrees) 


Fig.  2-The  scattering  function,  /( p),  as  deduced  from  the  sunlight 
observations  reported  in  Pyaskovskaya-Fesenkova  1957.  As  in  Figure  1, 
all  the  quantities  in  equation  (15)  are  known  so  that  the  scattering 
function  may  be  determined.  The  scattering  function  is  plotted  versus 
the  scattering  angle,  p,  for  sites  with  extinction  coefficients  of  0.15 
(crosses)  and  0.24  (dots).  The  scattering  function  from  equation  (21)  is 
drawn  as  the  smooth  curve.  Note  that  the  scatter  about  the  model  is 
small  and  that  there  is  no  significant  difference  between  the  two  sites,  as 
this  is  the  empirical  justification  that  the  functional  form  of  equation  (15) 
is  valid. 
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log-linear  plot.  From  this  plot  we  will  adopt  a formula  for 

/m(P)  = 108,15-p/4°  . (18) 

Now  the  total  scattering  function  can  be  found  as  a func- 
tion of  p from  equations  (16),  (17),  and  (18),  and  is  dis- 
played as  a curve  in  Figures  1 and  2. 

For  angular  separations  less  than  10  degrees,  the  data 
in  King  1971  show  that  equation  (18)  is  not  valid.  Schaefer 
1991  has  shown  that 

/m(P)  = 6.2  X 107  p-2  . (19) 

The  scattering  angle  is  measured  in  degrees  for  both 
equations  (18)  and  (19).  This  small  angle  scattering  forms 
the  aureole. 

The  reader  may  wish  to  compare  our  scattering  func- 
tion, using  appropriate  scaling,  to  that  of  McClatchey 
et  al.  1978  (Fig.  26  in  Section  14)  and  that  of  Rozenberg 
1966  (Fig.  72). 

5.  Discussion 

The  final  equation  for  Bmoon  is  derived  from  equa- 
tions (3),  (8),  (9),  and  (15)— (19).  For  Moon/sky  separa- 
tions greater  than  10  degrees,  these  results  can  be  sum- 
marized as 


j*  _ -^0  -0.4(3.84  + 0.026  lal +4  X 10-9  «4) 

(20) 

/( p) 

= 10536  [1.06  + cos2(p)l  + 10615-p/4°  , 

(21) 

X(Z)  = (1  - 0.96  sin2Z)-05  , 

(3) 

Braoo„=/(P)/*10-°'4W(Z”)[1-I0-0-4W<Z>]  • 

(15) 

The  required  input  for  this  general  formulation  is  the 
lunar  phase  angle  a,  the  Moon/sky  separation  p,  the 
extinction  coefficient  k , the  zenith  distance  of  the  Moon 
Zm,  and  the  zenith  distance  of  the  sky  position  Z.  These 
equations  are  not  as  formidable  as  they  might  look.  A 
pocket  calculator  can  evaluate  Bmoon  in  about  a minute  and 
only  a few  lines  of  code  are  needed  in  a computer  pro- 
gram. 

The  change  in  the  V-band  sky  brightness  caused  by  the 
moonlight  will  be 

AV  = -2.5  log  [(Bmoon  + B0(Z))/B0(Z)]  . (22) 

As  a guide,  we  have  tabulated  typical  values  of  Bmoon  (in 
nL)  and  AV  (in  magnitudes  per  square  arc  second)  for  a 
range  of  a and  p (see  Table  2).  This  table  was  constructed 
with  the  assumption  that  k = 0.172  mag/air  mass  (the 
median  V-band  extinction  at  the  2800-m  level  of  Mauna 
Kea),  and  that  Vzen  = 21.587  mag  sec-2  (the  mean  dark-sky 
value  at  the  same  site,  equivalent  to  Bzen  = 79.0  nL). 
Strictly  speaking,  the  values  of  Bmoon  and  AV  in  Table  2 
only  apply  to  a site  with  the  same  median  extinction 
coefficient  and  mean  dark-sky  brightness  as  the  2800-m 
level  of  Mauna  Kea  (e.g. , McDonald  Observatory). 

Given  the  two  extinction  terms  in  equation  (15),  a table 


TABLE  2 


Lunar  sky 

brightness 

effect  from 

model3 

Phase 

angle 

(a) 

5° 

Angular  distance  between 
moon  and  sky  position  (p) 
30°  60°  90° 

120° 

30° 

7216 
' (-4.48) 

1160 

(-2.87) 

530 

(-2.22)  ( 

437 

-1.93) 

818 

(-2.16) 

O 

O 

3364 

(-3.67) 

541 

(-2.13) 

247 

(-1.54)  ( 

204 

-1.30) 

381 

(-1.49) 

90° 

1351 

(-2.73) 

217 

(-1.35) 

99 

(-0.88)  ( 

82 

-0.71) 

153 

(-0.85) 

120° 

391 

(-1.58) 

63 

(-0.58) 

29 

(-0.34)  ( 

24 

-0.26) 

44 

(-0.32) 

aValues  of  Bmoon/  measured  in  nanoLamberts . For  lunar 
zenith  angle  of  60°  and  V-band  extinction  of  0.172  mag/ 
airmass.  For  this  table  p is  measured  along  the  great 
circle  on  the  sky  passing  through  the  moon  and  the  zenith. 
Therefore,  the  column  p = 60°  corresponds  to  the  zenith. 
Values  in  parentheses  are  AV  in  mag/sec2,  using  V = 

21.587  mag/sec2  (Bzen  =79.0  nL)  as  the  nominal  zenith  sky 
brightness,  and  scaling  the  zenith  sky  brightness  to  the 
nominal  value  at  the  zenith  angle  corresponding  to  p. 

comparable  to  Table  2 for  the  Mauna  Kea  summit  (with 
median  V-band  extinction  of  0. 113  mag/air  mass;  Krisciu- 
nas  1990)  gives  Bmoon  values  which  are  all  lower.  The 
extreme  example  would  be  for  a satellite  in  Earth  orbit, 
for  which  k = 0 and  Bmoon  always  equals  0 (i.e.,  the  sky  is 
black).  For  a dark  sky  site  with  k — 0.30  the  sky-bright- 
ness effect  due  to  the  Moon  is  larger  than  those  values 
tabulated  in  Table  2.  Because  the  moonlight  contribution 
comes  from  scattering  of  light,  where  one  is  observing 
through  more  atmosphere  (i.e.,  sea  level  vs.  on  a moun- 
taintop),  there  will  be  more  scattering. 

For  Palomar  and  Mount  Wilson,  which  have  V-band 
extinction  comparable  to  that  at  the  2800-m  level  of 
Mauna  Kea,  but  different  dark-sky  brightnesses  (due  to 
artificial  light),  one  can  use  the  values  of  Bmoon  from  Table 
2 and  equations  (1),  (2),  (3),  and  (22),  with  an  appropriate 
value  of  Vzen  for  equation  (1),  to  calculate  the  sky-bright- 
ness effect  (AV)  in  mag  sec-2. 

One  of  the  principal  uncertainties  in  using  equation 
(22)  on  any  given  night  comes  from  the  uncertainty  of  the 
adopted  dark-sky  zenith  brightness.  Even  if  it  is  mea- 
sured, one  does  not  know  the  exact  contribution  of  faint 
stars  in  the  beam  or  what  the  effect  is  of  the  solar  wind  on 
the  air  glow  on  that  specific  day.  See  Krisciunas  1990  for 
further  discussion. 

The  accuracy  of  our  general  formulation  can  be  deter- 
mined from  the  scatter  of  the  observed  brightnesses 
about  the  predictions.  In  Figure  3 we  plot  Bmoon  (ob- 
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Fig.  3-Observed  values  of  sky  brightness  due  to  scattered  moonlight 
(in  nanoLamberts)  vs.  model  values  from  equation  (15).  Open  circles: 
data  for  which  Bmoon  < B0(Z).  Squares:  data  of  1991  January  4.  Large 
dots:  all  other  data.  A line  of  slope  1 is  indicated, 
served)  vs.  Bmoon  (model)  for  the  data  of  Krisciunas  1990 
and  from  Table  1 of  this  paper.  If  we  take  the  ratios 

[Bmoon  (observed)  - Bmoon  (model) ]/Bmoor,  (observed) 

we  find  an  rms  variation  of  this  ratio  of  23%  for  the  Mauna 
Kea  data  (with  Bmoon  > B0(Z)),  8%  for  the  k = 0. 15  Russian 
site,  and  11%  for  the  k = 0.24  Russian  site.  Hence,  we 
conclude  that  our  formula  has  an  accuracy  of  8%  to  23%  in 
the  prediction  of  the  sky  brightness  from  moonlight.  The 
accuracy  will  be  poorer  when  the  Moon  or  sky  position  is 
close  to  the  horizon  so  that  small  uncertainties  in  the 
extinction  coefficient  will  be  magnified. 

To  the  best  of  our  knowledge,  no  formula  for  moonlight 
brightness  has  ever  appeared  in  the  literature.  Our  for- 
mula has  the  advantages  of  having  the  correct  functional 


dependencies,  yet  being  easy  to  use  and  accurate  to 
better  than  23%.  As  such,  the  equations  can  be  used  to 
optimize  exposure  times  for  photography  and  CCD  im- 
ages, to  predict  limiting  magnitudes  for  visual  observa- 
tions, as  well  as  to  predict  the  visibility  of  stars,  planets,  or 
comets  near  the  Sun  or  Moon  for  historical  purposes. 

We  thank  Ewen  Whitaker  for  useful  comments  and 
references  to  other  work. 
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